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Abstract 

Transmission of Dirac fermions through a chip of graphene under the effect of magnetic field 
and a time vibrating double barrier with frequency w is investigated. Quantum interference within 
the oscillating barrier has an important effect on quasi-particles tunneling. A combination of both 
a time dependent potential and a magnetic field generate physical states whose energy is double 
quantified by the pair of integers (n,l) with high degeneracy. The large number of modes that exist 
in the energy spectrum presents a colossal difficulty in numerical computations. Thus we were 
obliged to make a truncation and limit ourselves to the central (n = 0) and two adjacent side band 
(n = ±1). 
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1 Introduction 


Graphene [1] is a single layer of carbon atoms arranged into a planar honeycomb lattice. Since its 
experimental realization in 2004 [2] this system has attracted a considerable attention from both 
experimental and theoretical researchers . This is because of its unique and outstanding mechanical, 
electronic, optical, thermal and chemical properties [3]. Most of these marvelous properties are due 
to the apparently relativistic-like nature of its carriers, electrons behave as massless Dirac fermions in 
graphene systems. In fact starting from the original tight-binding Hamiltonian describing graphene 
it has been shown theoretically that the low-energy excitations of graphene appear to be massless 
chiral Dirac fermions. Thus, in the continuum limit one can analyze the crystal properties using 
the formalism of quantum electrodynamics in (2+l)-dimensions. This similarity between condensed 
matter physics and quantum electrodynamics (QED) provides the opportunity to probe many physical 
aspects proper to high energy physics phenomena in condensed matter systems. Thus, in this regard, 
graphene can be considered as a test-bed laboratory for high energy relativistic quantum phenomena. 

Quantum transport in periodically driven quantum systems is an important subject not only of 
academic value but also for device and optical applications. In particular quantum interference within 
an oscillating time-periodic electromagnetic field gives rise to additional sidebands at energies e + lhcj 
(/ = 0, ±1, • • •) in the transmission probability originating from the fact that electrons exchange energy 
quanta huo carried by photons of the oscillating field, cu being the frequency of the oscillating field. The 
standard model in this context is that of a time-modulated scalar potential in a finite region of space. 
It was studied earlier by Dayem and Martin [4] who provided the experimental evidence of photon 
assisted tunneling in experiments on superconducting films under microwave fields. Later on Tien 
and Gordon [5] provided the first theoretical explanation of these experimental observations. Further 
theoretical studies were performed later by many research groups, in particular Buttiker investigated 
the barrier traversal time of particles interacting with a time-oscillating barrier [6]. Wagner [7] gave 
a detailed treatment on photon-assisted tunneling through a strongly driven double barrier tunneling 
diode and studied the transmission probability of electrons traversing a quantum well subject to a 
harmonic driving force [8] where transmission side-bands have been predicted. Grossmann [9], on the 
other hand, investigated the tunneling through a double-well perturbed by a monochromatic driving 
force which gave rise to unexpected modifications in the tunneling phenomenon. 

In [10] the authors studied the chiral tunneling through a harmonically driven potential barrier in 
a graphene monolayer. Because the charge carriers in their system are massless they described the 
tunneling effect as the Klein tunneling with high anisotropy. For this, they determined the transmis¬ 
sion probabilities for the central band and sidebands in terms of the incident angle of the electron 
beam. Subsequently, they investigated the transmission probabilities for varying width, amplitude and 
frequency of the oscillating barrier. They conclude that the perfect transmission for normal incidence, 
which has been reported for a static barrier, persists for the oscillating barrier which is a manifestation 
of Klein tunneling in a time-harmonic potential. 

The growing experimental interest in studying optical properties of electron transport in graphene 
subject to strong laser fields [1 ] motivated the recent upsurge in theoretical study of the effect of 
time dependent periodic electromagnetic field on electron spectra. Recently it was shown that laser 
fields can affect the electron density of states and consequently the electron transport properties [12]. 
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Electron transport in graphene generated by laser irradiation was shown to result in subharmonic 
resonant enhancement [13]. The analogy between spectra of Dirac fermions in laser fields and the 
energy spectrum in graphene superlattice formed by static one dimensional periodic potential was 
recently performed [14]. In graphene systems resonant enhancement of both electron backscattering 
and currents across a scalar potential barrier of arbitrary space and time dependence was investigated 
in [ 5] and resonant sidebands in the transmission due to a time modulated potential was studied 
recently in graphene [16]. The fact that an applied oscillating field can result in an effective mass 
or equivalently a dynamic gap was confirmed in recent studies [17]. Adiabatic quantum pumping 
of a graphene devise with two oscillating electric barriers was considered [18]. A Josephson-like 
current was predicted for several time dependent scalar potential barriers placed upon a monolayer 
of graphene [19]. Stochastic resonance like phenomenon [20] was predicted for transport phenomena 
in disordered graphene nanojunctions [21]. Further study showed that noise-controlled effects can 
be induced due to the interplay between stochastic and relativistic dynamics of charge carriers in 
graphene [22]. 

Very recently, we have analyzed the energy spectrum of graphene sheet with a single barrier 
structure having a time periodic oscillating height in the presence of a magnetic field [23]. The 
corresponding transmission was studied as a function of the energy and the potential parameters. 
We have shown that quantum interference within the oscillating barrier has an important effect on 
quasiparticle tunneling. In particular the time-periodic electromagnetic field generates additional 
sidebands at energies e+lhuj(l = 0, ±1, • • •) in the transmission probability originating from the photon 
absorption or emission within the oscillating barrier. Due to numerical difficulties in truncating the 
resulting coupled channel equations we have limited ourselves to low quantum channels, i.e. I — 0, ±1. 

We extend our previous work [23] to consider monolayer graphene sheet in the presence of magnetic 
field but with double barriers along the x-direction while the carriers are free in the ^/-direction. The 
barrier height oscillates sinusoidally around an average value Vj with oscillation amplitude Uj and 
frequency u. The spectral solutions are obtained in the five regions forming our sheet as functions of 
different physical parameters. These are used to calculate the current density and therefore evaluate 
the transmission probability for the central band and close by sidebands as a function of the potential 
parameters and incident angle of the particles. We present our numerical results and discuss their 
implications for low quantum channels. 

The manuscript is organized as follows. In section 2, we present our theoretical model by defining 
the governing Hamiltonian and setting the applied potentials and external magnetic field. We solve 
the resulting eigenvalue equations to obtain the solutions of the energy spectrum for the five regions 
composing our system in section 3. Using the boundary conditions as well as the current density we 
exactly determine the transmission probability in section 4. Our main results and comparisons with 
existing literature will be presented in section 5. We conclude by summarizing our main results in the 
last section. 

2 Hamiltonian of the system 

Consider a two-dimensional system of Dirac fermions forming a sheet graphene. This sheet is subject 
to a vibrating double barrier potential in addition to a mass term and an externally applied magnetic 
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field as shown in Figure 1. Particles and antiparticles moving respectively in the positive and negative 
energy regions with the tangential component of the wave vector along the x-direction have translation 
invariance in the ^/-direction. Dirac fermions move through a monolayer graphene and scatter off a 
double barrier potential whose height is oscillating sinusoidally around Vj with amplitude Uj and 
frequency uj. The carriers are also subject to a magnetic field perpendicular to the graphene sheet 
B = B(x,y)e z and a mass term is added to a vector potential coupling. Dirac fermions with energy 
E are incident with an angle <j)\ with respect to the x-axis, the conservation of energy allows the 
appearance of an infinite number of modes with levels E + mhcj (m = 0, ±1, ±2 • • •). The Hamiltonian 
governing the system is composed by two independent terms (Hq,Hi) 

H = H 0 + H x (1) 


where the first part is 

Hq = vp(j • + -A(x, y)^j + V(x)l 2 + ft© (d\ — x 2 ) a z (2) 

and the oscillating barrier potential is defined in each scattering region by (see Figure 1) 

H{ = Uj cos (ut + SjO(x)). (3) 

The Hamiltonian H\ describes the harmonic time dependence of the barrier height, ft is the mass 
term, vp the Fermi velocity, a = (cr x ,a y ) are the usual Pauli matrices, phase difference 5j, I 2 the 
2 x2 unit matrix, the electrostatic potential V(x) — Vj in each scattering region and the magnetic 
field B(x,y) — B{x). Adopting the Landau gauge which allows the vector potential to be of the form 
A = (0 ,A y (x)) with d x A y (x ) = B(x), the transverse momentum p y = — id y is thus conserved. The 
magnetic field B = B$e z (with constant Bq) within the strip \x\ < d\ but B — 0 elsewhere 

B(x,y) = B 0 @(dl-x 2 ) (4) 


with the Heaviside step function © 


©(*) 


1, x > 0 

0 , otherwise. 


( 5 ) 


The static square potential barrier V(x) is defined by its constant value Vj in each region, similarly 
for the amplitude of the oscillating potential Uj 


V(x) = V :t = { 


V 2 , 

Vs, 


d\ < \x\ < G ?2 
\x\ < d\ 


Uj = < 


C/ 4 , d\ < x < d 


0 , otherwise 


0 , otherwise 


( 6 ) 


where the index j = 1, 2, 3,4, 5 denotes each scattering region as shown in Figure 1. 

Concerning the applied magnetic field, it is a constant and uniform magnetic field B perpendicular 
to the graphene sheet but confined to a strip of width 2d. Due to incommensurate effect and interaction 
with substrate, graphene can develop a mass term in the Hamiltonian. The vector potential that 
generates our magnetic field can be chosen of the following form 



—di, x < d\ 

x , \x\ < d\ (7) 

c?i, x > d\ 


with the magnetic length defined by Ip 



in the unit system {h — c — e— 1 ). 
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U4 cos(uit+6) 


V4 


Figure 1: (Color online) Schematic of a graphene monolayer in the presence of an oscillating potential 
and a magnetic field. Different scattering regions are indicated by an integer j=l,2,3,4,5. 


3 Spectral solutions 

We emphasize that the system Hamiltonian (1) is composed of two sub-Hamiltonian, H\ plays the 
role of a perturbation term with respect to Hq. The independence of these Hamiltonians leads to their 
commutativity [Hq, H\\ — 0 and therefore the corresponding eigenspinors ^ are the tensor product of 
two eigenspinor t/iq and 1 associated with Hq and respectively i.e. ^(x^y^t) = y)^(t) and 

the eigenvalue of H is the sum of eigenvalues E = Eq + E\. The eigenspinor ^ of the system obeys 
the equation 

H'ip(x,y,t) = id t ipj(x,y,t) (8) 

which can be written as 


[E 0 + Uj cos (cut + 5jQ(x))]ipj(x,y,t) = id t ipj(x,y,t). (9) 

The integration between to = 0 and t gives 

iPj(x, y, t ) = ipj( x i y, 0)e- iEot e~ iU i sin M+« j 0 (*))/« ( 10 ) 

where the last term is in the form of e iaj sm ®, which can be expanded into trigonometric series as 


+ OO 


gZ Oij sin <E> 


E 




(ii) 


with Oij — -jj- and = uot + 5jQ(x). Hence finally we obtain 

+oo 

^ j (x,y,t) = i; j (x,y, 0)e~ iEot E 

m =—oo 

and C m = J m (aj) satisfies the recurrence relation 

2l7lCm = &j (Cra+1 T Cm— l) 


( 12 ) 


( 13 ) 


where J m is the m-th order Bessel function of the first kind. Using these eigenspinors we readily 
determine the total energy from ( 12 ) to be 


E = Eq + mu. 


(14) 
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Taking into account energy conservation, the wave packet that describes our carrier in the j-th region 
can be expressed as a linear combination of wave functions at energies Eq + Iuj (l = 0, ±1, • • •). This 
is 

m,l =+oo 

y, t ) = e ikyV Y, $(*> y) J m-l («i) (15) 

m,Z=—oo 

where we have set e = ^ and w = Subsequently, the spinor (x, y) will be determined in each 
region j. 

The Dirac eigenvalue equation in the absence of oscillating potential for the spinor xp(x , y ) = 
('ip + ,'ip-) T at energy Eq reads 

HqiKx, y) = E 0 ip(x, y). (16) 


Using the explicit form of Hq given by (2) we find 


Vj + /i 

Pjx + ip y + iA(x) 


Pjx - ipy - iA{x) \ / VU 

vj- p J \ 8- 


= 6 


*/>+ 


(17) 


Vi 


where Vj = ^ and \± = ^ . Due to the translational invariance along the ^/-direction, the two- 
component pseudospinor can be written as rp±(x,y) — (p±(x)e lkyV . In region j — 1, 2, 4 and 5, we 
easily obtain the following two linear differential equations 


(Pjx -Py- iA(x)) if- = (e- Vj)cp + 
(Pjx+Py + iA(x)) <p+ = (e- Vj)(p-. 


(18) 

(19) 


In accordance with (15), the general solution in the j-th scattering region reads as 

1 


771 , 1 =+OO 

rpi (x,y,t) = e ikyV £ 


-1 

p 

o- to. 

li ) 


v z l J 


e ik i (as—ajj) + Jyj 


e-irfix-Xj) 


m,l =—oo _ 

XJm-l {ptj) e -^m-l)8 je -iv F {e+m^)t 


( 20 ) 


where sj = sgn (e + lm — Vj), the sign again refers to conduction and valence bands, Xj are the positions 
of the interfaces (Figure 1): x\ — X 2 = — g? 2 5 ^3 = — di, £4 = di, x$ = g? 2 - Note that, outside the 
barrier regions where the modulation amplitude is ay = 0 we have the function J m -l (ay) = The 
wave vector is given by 


which leads to the corresponding eigenvalues 


k l = s \\K e - Vj + lw ? - (ky + ^r 


e-Vj + lm = sj^(kj^J + (ky + jt 


( 21 ) 


( 22 ) 


with the magnetic length defined by Ib — yJl/B$ and the complex parameter 


IS 


7 3 - 
^ — S l 


k{ +i(ky + jfj 


V / ( t? ) 2+ ( 


U -u A 
1 /2 

L B 


= sje 1 ^ 


(23) 
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<$ = tan 1 (k y /kj) and the parameter d is defined by 

d = 


d\, x < —di 
—d\, x > d\. 


(24) 


Let us proceed to write down the solution in the intermediate zone j — 3 {—d\ < x < d\) contain¬ 
ing the mass term in addition to a perpendicular magnetic field. To diagonalize the corresponding 
Hamiltonian we introduce the usual boson operators 


ai - (<9 X . + ky + -p 


t Ib 
1 = ~F 

1 V2 


„ 7 X 

— 0 X + Gy + -7T- 


(25) 


which satisfy the commutation relation 


ai,a[ 


Sik • In terms of a/ and a|, equation (17) reads 


^3 + V ~ij§'<H 

l i} a \ ^3-M 


^ ' =(e + lw)l m 


‘Pi, 2 


‘-Pi,2 


(26) 


or in its explicit form 


= (e + lw -v 3 - n) ip t l 
Ib 

V2 t 

i—a] w = (e + Im - v 3 + n) 


Combining the above equations, we obtain for 


((e + Iw - v 3 ) 2 - /j , 2 ) ipi t i = J- aia\(pi tl . 


(27) 

(28) 

(29) 


It is clear that is an eigenstate of the number operator N = a\ai and therefore we identify 

(30) 


with the eigenstates of the harmonic oscillator \n — 1), namely 

<Pl ,l H n - 1) 


and the associated eigenvalues are 


e — v$ + lm = \j (/i/s) 2 + 2 n. 

I>B V 


(31) 


Finally, the solution in region j — 3 can be expressed in accordance with equation (15), as follows 


l —-|-oo 


M x ,y,t) = e ikvV E {aM + b!<PT)e- iMe+llByt 


(32) 


where cpf are given by 






± 


l = — C 


eJ t^ D [(e,, n l B )*-(»l B y}/2-l [±V^ (lt(* - ^3) + V*) 




(33) 


/ 


In the forthcoming analysis, we will see how the obtained results so far can be applied to deal with 
different issues. More precisely, we will focus on the transmission probability for different channels. 
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4 Transmission probability 


Based on different considerations, we study interesting features of our system in terms of the corre¬ 
sponding transmission probability. Before doing so, let us simplify our writing using the following 
shorthand notation 

At = , (34) 


B tn = 


y/%n(£l,n ± M) 

* 1,1 = D [{e,,nl B )*-UdB)*\/ 2 -l [±V^(Mb 

tl,l = D [(ei,nlB) 2 -(lAl B ) 2 }/‘2 ±V 2 (kyl B )- 

± /- ( 2d\ 

^ 2 ,Z [(o,n^) 2 -(^s) 2 ]/ 2-1 l 


^ 2 ,Z ^[(e^ n / jB ) 2 -(/x / B ) 2 ]/2 


i + kylB 

±V2 (l+v,)' 


Realizing that j e zrm; F^ j are orthogonal, we obtain set of simultaneous equations emanating from 
the boundary conditions at x = — g?2 

1=00 


+ - X! (°^ + tf) 

l=—oo 

1=00 / -. \ 

4 , -= E (“?*? - 6 ?ip) • 7 ™-< S) 

,n l— — r vn x t 7 


similarly at rr = — c?i 


l = OQ 

7 3 4+ 00+ -L ft 3 4- r7 - _ ^2 ikf(d 2 -di) , l2 p -ikf(d2-di)\ j , ( — \ 


l=—oo 
1=00 


1 = 00 / 1 \ 

3 ,n+ C+ _ j,3 D- t— _ \ ' | /7 2 r 2 p ikf(d,2-di) _ l2_ -ikf(d,2-di)\ t ( ^2 \ 

■m^m.nSl.m °m ll5 m,n l il,m ~ 2_> \^l Z l e °l z 2 e j^Kvo) 


^ (43) 


and at x = c?i 


1=00 

7 3 4 + 77+ I ft 3 4 - 77- — X~^ ^4 , l 4 \ t -z( m -Q (5 

l m z± m,n'l2,m ' u m^ ± m,n'l2,m / j \ a l ' u l ) J rn—l y y ^ 

l=—oo 

l —00 / 1 \ 

~ b m iB m,n^2,m = Y2 (“) e~ l{m ~ l)S . 

/ = — rxr> x l ' 


However, at x = o?2 we have 

1=00 

a m + + = X/ 

1 = —OO 


b b Lb _ 

a m z m u m 5 

Zrrn 


Jm—l ( 


| e -i(m-l)6 

(46) 


\W ) 


-ikf(d 2 - 

-di) N 

y 


(47) 
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As Dirac electrons pass through a region subject to time-harmonic potentials, transitions from 
the central band to sidebands (channels) at energies e ± mm {m = 0,1, 2, • • •) occur as electrons ex¬ 
change energy quanta with the oscillating field. It should be noted that (40-47) can be written in a 
compact form as 


Mil M12 
M21 M22 


^5 

^5 


^5 

^5 


(48) 


where the total transfer matrix M = M(l, 2) • M(2, 3) • M(3,4) • M(4, 5) and M(j. j + 1) are transfer 
matrices that couple the wave function in the j-th region to that in the (j + l)-th one. These are 
explicitly defined by 


-l 


M(l, 2) = 


M(2, 3) = 


M(3,4) = 


M(4, 5) = 


N+ N" 


¥+ 

y 2 - 


Ft Fn 


-1 


-1 


Kf K: 


-1 


K+ 


c 

C 

G+ 

G“ 

Qt 

Qr 

F^ 

F r 

Di 

Di 

D+ 

D 2 

I 

I 

E+ 

E“ 


whose matrix elements are expressed as 

( C W = J m-l (£) 

( G± )m,Z = Hzi) ±l Jm-l (^) 

( Y f)w =e±ifc?(d2_dl)j ^ (|) 

(*2 ) m> , = ±(zf) ±1 e ±ik i^~ d ^J m _ l 0 ) 

(Qr ) m: i = A m,nVt,m 8 m,l 
( F r ) m j = 


(» l)m,l = Jm-l g) e ( - i(m “ 05) 


4' 


Z )m,l = (^) 

= e ±ik t^~ dA J m _, 0) 

= ±(^) ± 1 e ±ifc ^ d2 - dl ) J m _, 0 ) 

( E 4 m i = ±(4) ±:l <W 


(49) 


(50) 


(51) 


(52) 


(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 
(62) 

(63) 

(64) 


and the unit matrix is denoted by I. We assume an electron propagating from left to right with 
quasienergy e. Then, r E {1,2}, Si = {a^} = {5 m ,o} and S 5 = { b is the null vector, whereas 
S 5 — {&m} — {tm} and S{ = {bjn} = {r m } are vectors associated with transmitted waves and 
reflected waves, respectively. From the above considerations, one can easily obtain the relation 




S 5 = (Mil) .=1 

which is equivalent to the explicit form 

/ t~N \ 


t -1 
£0 
h 


\ f N ) 


= (Mil) 


-1 


/ 0 \ 
0 
0 
0 
1 
0 
0 
0 

V° J 


(65) 


( 66 ) 


The minimum number N of sidebands that need to be taken is determined by the strength of the 
oscillating potential, N > max(^,^) [10]. Then the infinite series for the transmission T can 
be truncated considering only a finite number of terms starting from — N up to TV. Furthermore, 
analytical results are obtained if we pick up small values of 0 L 2 — oq = ^ and include only the 
first two sidebands at energies e ± mm along with the central band at energy e. This gives 


t-N+k — M [k + 1, N + 1)] 


(67) 


where k — 0,1, 2, • • • , 27V and M.' denotes the inverse matrix (Mil) -1 . 

Using the reflected J re f and transmitted J tra currents to write the reflection and transmission 
coefficients Ri and Ti as 


Ti 



Ri 


| ^inc,0 


( 68 ) 


where Ti is the transmission coefficient describing the scattering of an electron with incident quasienergy 
e in the region 1 into the sideband with quasienergy e + lw in the region 5. Thus, the rank of the trans¬ 
fer matrix (M) increases with the amplitude of the time-oscillating potential. The total transmission 
coefficient for quasienergy e is 

Z=+00 

T = T t . (69) 


l = — OO 

The electrical current density J corresponding to our system can be derived to be 


which explicitly reads as 


J = v F ^a x ip 


J\ n c,0 = VF (4 + i z oY) 

Jref.l = V F (b[ ) * b) (z/ + (z, 1 )*) 

Jtra,l = V F (of)* of (zf + (zf)*) . 


(70) 

(71) 

(72) 

(73) 
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The transmission coefficient for the sideband, TJ, is real and corresponds to propagating waves. It can 
be written as 


Ti = 


sfkf 




{k o ) 2 + (k y ~ fr'j 
(kf) 2 + (k y + f|) 


j.5 |2 


(74) 


Now using the energy conservation to simplify 7) to 


Ti 


k± 

k o 



j .5 |2 

I • 


(75) 


To explore the above results and go deeply in order to underline our system behavior, we will 
pass to the numerical analysis. For this, we will focus only on few channels and choose different 
configurations of the physical parameters. 


5 Discussions 

We discuss the numerical results for both the reflection and transmission coefficients, which are shown 
in Figures 2, 3, 4, 5, 6 , 7, 8 , 9 for different values of the parameters (e, v 2 , ^ 3 , di, g? 2 5 04 , 5, /x). 

To start with we point out the efficiency and accuracy of our computational method and compare 
our results with those reported in the literature. As a matter of fact, Figure 2(a) and Figure 2(b) 
reproduce exactly the results obtained in [24] for single barrier and [25] for double barrier, respectively, 
with the proper choice of parameters. Note that reference [24] was the first to treat the confinement 
of Dirac fermions by an inhomogeneous magnetic field. These polar graphs show the transmission as 
a function of the incidence angle, the outermost circle corresponds to full transmission, Tq = 1 , while 
the origin of this plot represents zero transmission, i.e. total reflection. For energies satisfying the 
condition (e + Ivj)Ib < we obtain total reflection [ 4,25]. This is equivalent to the condition on 
the incidence angle (j) < (j) c where </> c is the critical angle given by 

^ <76) 

which is analogous to the case of light propagation from a refringent medium to a less refringent one. 

After a satisfactory confirmation that our numerical approach reproduced published results, we 
plot the transmission versus the phase shift S in the presence temporal barrier oscillations for the cases 
a — 0 L 2 — 04 and a = 0 L 2 — Figure 3 illustrates the variation of the transmission coefficient To as 
a function of the phase shift S for various harmonic amplitudes a (0 < a < 0.99). The first point to 
emphasize is that both plots in Figure 3 are periodic with period 2 tt, which is obvious. For a = 0, the 
transmission does not depend on the phase shift since the vibration amplitude has been set to zero. 
For a = 0 L 2 — &4, the transmission varies sinusoidally between 0 and 1, the maximum for different 
plots does not change for various values of a (Figure 3(a)). However, when the harmonic oscillation 
amplitudes are not equal (04 / 04 ), we observe that there is a remarkable change in the evolution 
of Tq versus S. We illustrate this situation by selecting 04 = 204 then one can see that as long as 
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Figure 2 : (Color online) Polar plot showing transmission probability (transmission Tq (Z=0)) as a function 
of angle </>J for different values of the parameters, (a): = {0.5,1.5, 3, 3.67}, els = 3.7, OL 2 — «4 = 0, 

6 = 0, V21b = vaIb = 0, fils = 0 and ^ 3 ^ = 0. (b): els — {0.6, 2,4, 8 }, ^ 0.5, = <^4 = 0, 

5 = 0, = 0.6, v^Ib — V 4 Ib = 0.5, fils = 1 and ^3/^ = 0.4 . 

such difference increases we observe a drastic change in the transmission from full transmission (total 
transmission) to zero transmission (total reflection), see Figure 3(b). 

T o T 0 




Figure 3: (Color online) The transmission coefficient Tq as function of phase shift 5 through graphene 
double barriers for fixed values cIb, kyls, I^Ib, ^ 3 Ib, b and V 2 Ib = ^ 4 Ib but for different 

values of a. We used els — 25, k y ls = 2, ^ = 0.3, — 1.35, /iIb — 4, v^Ib 4, vjIb — 2 and 

V 2 Ib = ^4 ^b = 6 and a varies from 0 to 0.99. (a): a — 0 L 2 — a 4 , (b): a — 0 L 2 — 2 <24 

Let us now demonstrate through Figure 4 how the first sideband transmissions T\ and T_i vary as 
function of the phase shift S. It is clearly seen that the central transmission To behaves sinusoidally but 
at some value of a we observe that Tq changes its behavior and becomes sharply peaked. T± and T_i 
show also sinusoidal behaviors with non-symmetric double humps in regions where Tq is suppressed. 
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However, one can see that there is a symmetry between the two double humps of the transmissions 
T\ and T—\. 


To , 1,-1 



Figure 4: (Color online) Graphs depicting the transmission probabilities as function of phase 5 for graphene 
double barriers with wIb — 8, ot 2 — «4 = 0.8, els 25, k y ls = 2, ^ = 0.3, ^ = 1.35, fils = 4, i^Z# = 4 
and V2Ib — ^4 Ib — 6. 



Figure 5: (Color online) Graphs depicting the transmission probabilities as function of potential v^Ib for 
the monolayer graphene barriers with = 1-5, = 0.5, V 2 Ib = ^ 4 Ib = 4, cIb — 25, k y ls = 2, /iIb = 4, 

hjZ# = 2, C42 = 04 = 0.5 and S = {0, |,7r, ^}. To (color blue), T_i (color red) and Ti (color green). 

Now we will study how the three bands: central and two lateral ones, vary depending on the phase 
difference of the oscillating potentials in the intermediate region (Figure 5). For different phase shifts, 
the transmissions of side bands T\ and T_i are dominant either before or after the bowl centered 
region in the propagation energy eZ#. In the vicinity of this bowl, one of the two transmissions is more 
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symmetrical with respect to a vertical axis passing through the energy cIb • The degree of dominance 
of the transmissions T\ and T_i is less pronounced in the case of advanced phase quadrature S = | 
(Figure 5(b)). But they are more dominant in the case of 5 = 0 and 5 = ^ where the potential v^Ib is 
greater than the propagation energy cIb and are less dominant when v^Ib is less than the propagation 
energy e/#. The behavior of the transmission side bands differs in the case of opposite phase shift 
S = n. On the other hand, the transmission of the central band, for different phase shifts, has also a 
dominance of either side of the high potential v^Ib than the propagation energy els or the other side 
where the potential v^Ib is small than the same energy. 



Figure 6 : (Color online) Graphs depicting the transmission probabilities as function of potential V 2 Ib for 
monolayer graphene barriers with ol 2 — &4 = 0.5, zuIb 2, els 25, k y ls = 2, = 0.5, = 1.6, 

tils = 4, ^ 3 Ib = 4 and 5 = {0, §, 7 r, ^}. To (color blue), T_i (color red) and T\ (color green). 

Figure 6 is similar to Figure 5, the only main difference to be noted is that there is presence of 
peaks in the bowl centered around the value V 2 — e and To ? _i 5 i heights are in the same order either 
before or after the bowl. These peaks are due to the resonances between the bound states existing 
in both sides of the regions subject to the potential ^ 3 . This behavior is normal if we keep in mind 
that our double barrier is composed of two successive squares with the same potential V 2 and width 
(g ?2 — d\) separated by the width 2 d\ corresponding to the intermediate region. 

Figure 7 shows the transmission probability as a function of the incident energy of electrons for 
v> 2 Ib = 25, kyls = 2, ^ = 0.5, — 1.5, fils = 4, v^Ib — 4 and different amplitudes of the oscillating 

barrier without shift 5 = 0. Resonant peaks are narrow and could have important applications in 
high-speed devices based on graphene as has been suggested previously [ L0] . The evolution of the 
central and two lateral transmission bands depend on the width of the double barrier potential over 
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time accompanied by a magnetic field (Figures 7(a), 7(b)). 



Figure 7: (Color online) Graphs depicting the transmission probabilities as function of energy els for 
the monolayer graphene barriers with 0 L 2 — &4 = {0.08,0.5}, V 2 Ib = 25, k y ls = 2, = 0.5, — 1.5, 

fils = 4, ^ 3 Ib = 4, tu/# = 2 and 5 = 0. 


7o.l-i 


F 0 i- 
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eh 


Figure 8: (Color online) Graphs depicting the transmission probabilities as function of energy els for 
monolayer graphene barriers with <^2 = ^4 = 0.99, = {0,0.02,0.5,1,1.3}, V 2 Ib = 25, kyls = 2, 

p- = 1.5, ills = 4, ^ 3 Ib = 4, ccjZ# = 2 and 5 — n. Tq (color blue), T_i (color red) and Ti (color green). 


Figure 8 presents transmission versus the system energy for different widths. Indeed, we observe 
that in Figure 8(a) as long as the width is very small the central band is dominant and therefore 
the transmission becomes total independently of the applied potential. Figure 8(b) is obtained by 
increasing the width d up to some value, one can see the dominance of the two sideband transmissions 
compared to central band one. We notice that these two sideband transmissions are symmetrical with 
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respect to an axis of symmetry located at double barrier potential V 2 of the propagation energy. After 
increasing the width d, we end up with Figure 8 (c), which is similar to the last one but this time 
with dominance of the central band transmission. It is clearly seen that the total transmission is 
less than or equal to unity. In Figure 8 (d), the central band transmission recovers its dominance but 
evanescence of two sideband transmissions. 



Figure 9: (Color online) Graphs depicting the transmission total T as function of energy gap / iIb for the 
monolayer graphene barriers with: ^ = 0.7, v^Ib — v^Ib — 4, cIb = 25, kyls = 2, v^Ib — 4, voIb — 2, 
a 2 = a 4 = 0 and 5 = 0 (a): g- = {0.09,0.2, 0.5}, a 2 = a 4 = 0 and S = 0. (b): g = 0.5, a 2 = a 4 = 0.7 
and 5 = ^ . 


Figure 9 is intended to see the influence of increasing the width of the intermediate zone, where 
there is a magnetic field, on the dominant transmission central depending on the mass term /. lIb 
that in the intermediate region. The distance c ?2 remains constant which means that the widths of 
regions 2 and 4 decrease if d\ increases. Figure 9(a) shows that progressively as the distance d\ 
increases, the central transmission acquires resonances which clamp by increasing amplitudes whose 
upper peaks correspond to a total transmission (maximum). The maximum value of To(/i) is the 
unit since a = OL2 — oq = 0. In Figure 9(b) for a 7 ^ 0 the maximum value of To(/i) decreases at 
the expense of transmission sidebands T_i ? i(/i). We note that the sum of the three transmissions 
T 0 ,_i,i(/i) converges whenever towards unity. 


6 Conclusion 

In this present work, we studied the transmission probability in graphene through double barriers 
with periodic potential in time. The double barrier contains an intermediate region has a magnetic 
field with a mass term, but the two temporal harmonic potentials with different amplitudes and phase 
shifted are applied one hand and on the other in both regions restricting the intermediate region. 
This panoply of potential makes our studied system rich in terms of physical states whose energy is 
doubling quantified by the pair (n, l) extensively degenerated with a very large number of modes. 

To identify the difficulties posed, we made the problem by adequate truncation to reduce all modes 
in three modes one central and two lateral indexed by (0, —1,1). We tried to study the influence of 
various parameters such as (e, V 2 , ^ 3 , di, k) on the transmission probability and 

highlight some properties of the system under consideration. The critical angle (see (76)), at which 


15 


















total reflection sets in, showed in an efficient manner the analogy between the propagation of Dirac 
fermions in our system and the propagation of the light of a more refractive homogeneous isotropic 
transparent medium a less refractive one. This built an interesting bridge between two areas of physics 
such matter and light. 

The transmission probability Tq(5) is obtained to be harmonic with frequency proportional to 
w of time dependent amplitude a. We observed that as long as the amplitude a of time-harmonic 
potentials is increased Tq(/i) is decreased at the expense of lateral transmissions T_yi(/i) and the three 
transmission behaves in a complementary manner and are bounded. While the sum of the 

three transmission converge towards unity, as required by the unitarity condition. 
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